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Abstract 

In an effort towards understanding turbulent scalar mixing, we 
study the effect of molecular mixing, first in isolation and then ac- 
counting for the effects of the velocity field. The chief motivation for 
this approach stems from the strong resemblance of the scalar prob- 
ability density function (pdf) obtained from the scalar field evolving 
from the heat conduction equation and that evolving in a turbulent 
velocity field. However, the evolution of the scalar dissipation is 
different for the two cases. We attempt to account for these differ- 
ences, which are due to the velocity field, using a Lagrangian frame 
analysis. After establishing the usefulness of this approach, we use 
the heat-conduction simulations (HCS), in lieu of the more expen- 
sive direct numerical simulations (DNS), to study many of the less 
understood aspects of turbulent mixing. Comparison between the 
HCS data and available models are made whenever possible. It is 
established that the (3 pdf characterizes, quite well, the evolution of 
the scalar pdf during mixing from all types of non-premixed initial 
conditions. 


1 Introduction 


Diffusion of smoke from chimneys into the atmosphere and mixing of fuel 
and air in an aircraft engine are but two of numerous instances in which 
scalars are convected and diffused in a turbulent flow field. Environmental- 
ists and engineers alike are interested in understanding this mixing process. 
For many applications, a knowledge of the basic statistics of mixing suf- 
fices. However, when scalar mixing is accompanied by chemical reaction, 
a detailed understanding of the probability density function (pdf) of the 
scalars is required (Pope 1985). The direct numerical simulations (DNS) of 
Eswaran and Pope (1988) (referred to as EP henceforth) gave us the first 
close look at the details of the evolution of the scalar pdf in constant density 
isotropic turbulence. Since then, there have been other DNS studies of pas- 
sive scalar mixing and mixing with chemical reactions (Givi 1989). These 
DNS data, apart from enhancing our understanding of turbulent mixing, 
also provide a test bed for evaluating mixing models. The mapping closure 
model (Chen, Chen and Kraichnan 1989, Pope 1991, Gao 1991) has been 
strongly validated by DNS data and has emerged as the leading turbulent- 
mixing model. The /3-pdf model for scalar mixing (Girimaji 1991a, 1992) 
represents an instance of the direct use of DNS observations to formulate 
simplified closure models for problems of practical importance such as tur- 
bulent combustion (Girimaji 1991b). Despite these recent long strides, our 
understanding of turbulent mixing remains incomplete. 

A popular technique often used in the past to understand turbulent 
processes is the study of the process in related simpler systems (e.g., use 
of Burger’s equation rather than Navier-Stokes). In the same vein, we first 
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study the evolution of initially-random scalar fields evolving according to 
the heat conduction equation rather than the full scalar evolution equation. 
Then we account for the effect of the velocity field using a Lagrangian frame 
analysis where these effects are related to the material-element deformation 
characteristics of the turbulent velocity field. 


1.1 Motivation 


We study the evolution of a scalar field </>(x, t ) subject to the heat conduc- 
tion equation, 

d£ _ D ay 

dt dxidxi ’ 

rather than the full evolution equation, 


d£, gg _ D 

dt ' dxi dxidxi ’ 


( 2 ) 


where D is the coefficient of Fickian diffusion and u(x, t) is the velocity 
field. There are several a priori motivations to perform such a study to 
better our understanding of turbulent mixing. For example, irrespective of 
whether the scalar field evolves according to equation (1) or (2), the scalar- 
pdf (F(<p)) evolution equation is the same. In either case, the pdf evolves 
in an isotropic turbulent field according to (Pope 1985) 

<9F(V>) d 2 


dt 


dt/>dip 


{FWxi'l’)}, 


( 3 ) 


where ?/> is the probability-space value of the scalar concentration and the 
conditional scalar dissipation xii*) i s given by 


XW = D{ 


dx t dxi 


\<i> = i>)- 


( 4 ) 
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The notation (b\r) is used to denote the conditional expectation of b with 
respect to r. For the purpose of analyzing the evolution of the scalar pdf, 
it is clear from equation (3) that time t can be simply rescaled to account 
for the diffusion coefficient and the diffusivity can be assumed to be unity 
without loss of generality. 

A strong a posteriori motivation for this study comes from our obser- 
vation that the pdf’s calculated from equation (1) (for the case of equal- 
quantity, non-premixed initial condition) at various stages of scalar decay 
are quite similar to those obtained from the DNS study of EP at the cor- 
responding values of scalar variance. It should be recalled here that one 
of the principal findings of the DNS study is that the shapes adopted by 
the scalar pdf are independent of the initial scalar length scale. That is, at 
any stage of scalar decay characterized by the value of variance, the shape 
of the scalar pdf is independent of the initial scalar length scale and tur- 
bulence statistics. At any given value of scalar variance, we find that the 
scalar pdf obtained from the heat conduction equation is quite similar to 
the scale-independent pdf in the turbulent field obtained by EP. 

Role of the velocity field. As a prelude to our study, we first attempt 
to understand the role of the velocity field in turbulent scalar mixing. The 
role of the velocity field can be most easily understood from the spectral 
version of equation (2), 

f rniiii(k — m,t)^(m,t)dm — k 2 ^(k,t), (5) 

ut J -oo 

where a tilde denotes the Fourier transform of the quantity and k is the 
wavenumber vector of magnitude k. From this equation we can glean that, 
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when the length scale of the velocity field (/„) is much larger than that 
of the scalar field ( l a ), the effect of the velocity field is negligible. On the 
other hand, when the length scales are comparable, or when the scalar 
length scale is larger, the velocity field cascades the scalar energy down to 
smaller scales until molecular diffusion ultimately takes over. However, this 
rearrangement of scalar spectrum is such that the scalar statistics in the 
physical space are unaffected. So, it is clear that in the former case when 
l u /l s > 1, the neglect of the velocity field is quite reasonable, whereas in 
the latter case (I u /l s < 1), the consequences of its omission are unclear. EP 
find that over a wide range of the initial scalar length scales (/„//„ ~ 1 to 
16) the evolution of the shape of the pdf is similar although the time scales 
of evolution are vastly different. This clearly implies that the velocity field 
has little influence on the scalar pdf shapes. Therefore, our present study 
of turbulent mixing using the heat-conduction equation (which represents 
l u /l 3 ~ * oo) should be meaningful at least in understanding the evolution 
of the shape of the scalar pdf. 

In the presence of a velocity field, the scalar gradient evolves according 
to 

d d<f> d d<f> _ dui d(j> d 2 d<j> . 

dt dxj U ' dxi dxj dxj dxi ^ dxfix i dxj 

The effect of the velocity field on the scalar gradient is not only to rearrange 

its spectral content, but also to affect its physical-space pdf, owing to the 

presence of the first term on the right-hand side (RHS). The evolution of 
the scalar-gradient pdf, F(g), where g is the probability-space counterpart 
of the scalar gradient vector, is given by 


= _ JL lF ( g)( _*!ii* + g jiig)i 

dt dgi ^ ' dxj dxi dx l dx x dxj 


(7) 
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So, it is clear that the velocity gradients directly affect the evolution of the 
scalar-gradient pdf and, hence, that of the scalar dissipation, 0 (= §£-§£-)• 
However, the effect of the velocity field on the scalar pdf is only indirect, 
via the scalar gradients. 

The scalar variance a 2 evolves according to 

da 2 d<j> d<j> 

~dt = “fete/ = ~ € ‘ 

where e s is the mean scalar dissipation. So, the velocity field affects the 
evolution of the variance via its effect on e a . 

1.2 Numerical Aspects 

The solution for equation (1) has long been known in integral form for gen- 
eral initial conditions and in closed form for many simple initial conditions. 
The object of this study is to perform detailed statistical analysis of the 
behavior of the solution of equation (1) evolving from random initial condi- 
tions. As mentioned before, we expect such a study to shed light on many 
aspects of turbulent mixing. 

In this paper, we restrict our attention to constant-density, isotropic, 
two-scalar mixing. Without loss of generality, we can assume the concen- 
tration of one of the scalars to be <f>(x,t) (0 < <f> < 1) and that of the 
other to be (1 — <f>). Since only one of the concentrations is an indepen- 
dent variable, it suffices to study the statistics of either one of the fields. 
In constant-density turbulence, scalar concentration and scalar mass frac- 
tion are equivalent and, hence, throughout this paper we use the terms 
synonymously. 
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The computational domain is a cubical box with periodic boundary 
conditions. The initial scalar field is random in nature. Owing to the 
complexity of the initial conditions, only a numerical solution to the heat 
equation is possible. The Fourier coefficient of each wavenumber evolves 
independently, and its value at any time t is given by 

4>{k,t) = 4>{ k,0)e-*\ (9) 

Hence, even time stepping is not necessary. Knowing the initial field, the 
Fourier coefficients at any required time t are determined using the above 
equation. Then, the scalar information is transformed to physical space and 
all the required statistics are evaluated. 

The initial conditions are specified in a manner similar to that of EP. 
The procedure involves three steps. 

1. First, the Fourier coefficients of various wavenumbers are assigned 
values according to 

<£*(k) = f(k)exp[2iria(k)], (10) 

where a (k) is a random variable distributed uniformly between zero 
and unity. The function /( k) is given by 

1 forfc s -3< k<k a + 3 (11) 

m = 

0 otherwise. 

The central wavenumber, k s , of the scalar energy containing range is 
determined from the prescribed initial scalar length scale according 
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to 



( 12 ) 


2. The scalar field in transformed into the physical space, yielding ^*(x). 

3. The initial scalar field <f>(x,t) is now assigned from the field ^*(x), 
depending upon the initial pdf required. For example, if we need the 
initial pdf to be 


0 ) = ijftW + <^(1 - 0 )], 

then the scalar values are assigned according to 


(13) 


1 if ^*(x) > 0, (14) 

<?H x , o) = 

0 if 4>*{x) < 0. 


A triple-delta initial condition can be obtained as follows: 


1 if 4>*{x) > q, (15) 

<^(x, 0) = 0.7 if ^ > <^*(x) > r, 

0 otherwise, 

which produces the initial pdf 

F(i>) = a x 6(ip) + a 2 S(ip - 0.7) + (1 - <n - a 2 )S(rp - 1), (16) 

where a x and a 2 depend upon the values of q and r. Thus, we can 
obtain an isotropic, random initial scalar field of prescribed length 
scale and pdf. 
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The number of scalar energy containing wavenumber nodes is different for 
different values of k 3 : the larger the initial scalar length scale, the smaller 
is the number of energy containing wavenumber nodes. Since there is no 
cascading mechanism in the heat conduction equation, the number of energy 
containing nodes is the same at all times for any given initial length scale. 

All the results presented in this paper are from simulations using either 
64 3 or 128 3 grid points. For the calculation of the scalar moments, it is 
found that 64 3 simulations are adequate. For the calculations of the pdf s 
and scalar dissipation statistics, 128 3 simulations are used. For the sake 
of accuracy, scalar gradients are first formed in Fourier space, then trans- 
formed to physical space, and finally used to calculate the scalar dissipation 
statistics. The scalar pdf’s were calculated from 200 equal-sized bins in the 
range (0, 1). 

Outline. First, we establish the similarities and the differences between 
molecular mixing with and without the turbulent velocity field. We achieve 
this by comparing the heat-equation data with the DNS data for the case of 
equal-quantity, non-premixed, two-scalar mixing. Then, we use these find- 
ings to effect simplifications in a Lagrangian frame analysis of turbulent 
scalar mixing - in an effort to account for the role of the velocity field. Af- 
ter establishing the usefulness of this approach, we use the heat-conduction 
equation to study the statistical behavior of scalar fields evolving from a 
variety of other important initial conditions. The data from these heat- 
conduction simulations (HCS) are then used to evaluate aspects of turbulent 
mixing models which have hitherto been untested, owing to the lack of cor- 
responding DNS data. Another objective of the paper is to determine if the 
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/? pdf, which was found to characterize the equal-quantity, non-premixed, 
scalar pdf quite well (Girimaji 1991a), can adequately replicate the scalar 
pdf for non-equal-quantity mixing also. 

The remainder of the paper is organized as follows. In Section 2, we 
study the equal-quantity, non-premixed mixing. In Section 3, we perform 
the Lagrangian frame analysis to account for the effects of the velocity field. 
In Section 4, we study non-equal (// ^ 0.5), non-premixed initial fields and 
partially-mixed initial fields. We conclude with a summary in Section 5. 

2 Equal- quantity, non-premixed mixing 

In this Section, we present the results of equal-quantity, non-premixed scalar 
fields evolving by the heat conduction equation. As the pdf evolves from 
the double-delta initial condition, the mean scalar mass-fraction (/i) remains 
constant at 0.5. Comparisons with the DNS of Eswaran and Pope (1988), 
the /?-pdf model of Girimaji (1991a), the scalar-dissipation model of O’Brien 
and Jiang (1991), and the scalar-dissipation model of Girimaji (1992) are 
made wherever appropriate. 

2.1 Evolution of scalar statistics 

In Figure 1, the pdf of scalar mass fraction obtained from the heat con- 
duction simulation (HCS, k s = 8), the DNS ( k a = 8) of Eswaran and Pope 
(1988), and the /9-pdf model of Girimaji (1991a) are compared at various 
stages of mixing. While comparing the results, it is important to recognize 
that the initial spectral fields of the HCS and DNS are slightly different. 
The high wavenumber components of the DNS initial field are modified to 
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yield somewhat smooth initial fields. The HCS initial field is unsmoothed. 
So, the DNS data at each level of scalar mixing may represent the pdf at a 
slightly more advanced stage of decay than the corresponding HCS or the 
model pdf. The jaggedness in the early-time HCS pdf is due to the presence 
of the untruncated high wavenumber components of the scalar field. Given 
all these factors, the agreement is quite good. Both HCS and the model 
capture the various phases of the DNS evolution well. 

The main inference from the DNS study is that, irrespective of the initial 
scalar length scale ( l a ), the shape of the scalar pdf at any given stage of 
decay is the same. This means that in turbulent mixing the normalized 
scalar moments are unique functions of the scalar variance, independent of 
l s . This brings up the question, are the normalized scalar- moments obtained 
from HCS also unique functions of sc ala r variance? In Figure 2a, the HCS 
flatness-factor is plotted as a function of variance (<r 2 ), for various k a ranging 
from 1 to 16. Except for the cases of k s = 1 and k s - 2, the plots do 

collapse to within statistical error. The flatness-factor of k a — 1 and k a — 2 
do collapse at the early stages of mixing but at later times show a small 
but distinct deviation to lesser values. This is in disagreement with the 
DNS results which indicate that the pdf shapes are the same for the entire 
range of k a (1 through 16). The difference can be attributed to the effect 
of the velocity field. As was seen in the previous Section, the scalar-energy 
cascade (due to velocity gradients) is particularly important for small k a . 
In real turbulence, the scalar spectrum for these cases would shift to larger 
wavenumbers, and the evolution subsequently would be the same as that of 
larger k a . As for larger k a , the flatness-factor asymptotes to a Gaussian value 
of 3 during the latter stages of mixing. It remains nearly 3 for about three 
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decades of variance and then falls off to smaller values for cr 2 > 10 -5 (not 
shown). This behavior was also observed by Gao (private communication) 
in the DNS of turbulent scalar mixing. Furthermore, at very small variance 
values, the statistical and numerical errors involved become too large to 
compute the high-order moments accurately. In any event, this completely- 
mixed limit is of lesser practical value and, hence, has not been investigated 
further. In Figure 2b, the super-skewness of the same cases as in Figure 2a 
are compared. Again for k a > 2, we have collapse. During the latter stages 
of mixing (10~ 5 < cr 2 < 10 -3 ), the super-skewness is close to the Gaussian 
value of 15. 

The main inferences from Figures 1 and 2 are 

1. The HCS captures the details of the DNS scalar pdf evolution for 
k a =8 (Figure 1). 

2. For large enough k a (> 2), the HCS-pdf evolution is universal (i.e., 
same for all k s ; see Figure 2). The universal evolutions of HCS and 
DNS must be the same, considering the good agreement between HCS 
and DNS for the case of k a = 8. Hence, the velocity field does not 
affect the shape of the scalar pdf for large enough k a (small l,). 

3. For k a < 2, the velocity field affects the shape of the scalar pdf but 
modifies it only slightly. For most practical purposes, this modifica- 
tion can perhaps be neglected. 

4. The /?- pdf model captures the universal evolution remarkably well. 
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Evolution of <j> m i n and <f> m ax> The quantities <j> m i n and 4> ma x are the 
global minimum and maximum values of the scalar concentration at any 
given time. The evolution of these are of particular interest in the context 
of the evolution of conditional scalar dissipation. It is shown in Girimaji 
(1992) that the values of conditional dissipation at ^ min and <f> max are zero. 
The pdf value is zero outside the range (<j> m in, <f>max)- The conditional scalar 
dissipation model of O’Brien and Jiang (1991) implies that <^ mtn and <f> m ax 
do not change in time. 

For the problem under consideration (periodic scalar field with initial 
double-delta pdf), the exact solution is 

-I 1 ±oo 

<f>(x,t) = - + —j= ^(k,0)exp[-fc 2 t]exp[v /= Tk.x]dk, (17) 
2 V27T- !*(-! 

where <^>(k, 0) is determined from the initial conditions. At time t 0, the 
(j) min and <f>max are clearly zero and unity. With time, the value of <j) at all x 
locations changes. The fastest change comes about at a location which has 
the highest small-scale contribution. The slowest change is at a location 
where the large-scale contribution is the highest. From this we can surmise 
the following: 


<t>max{t) < £ + ^exp{~t], 

(18) 

Km{t) > | - xp[-i ]. 

(19) 


In Figure 3, the temporal evolution of ^> m ,„ and <f> ma x for various k s are 
plotted. The rate of decay of <f> mar increases with k„. Similarly, the growth 
rate of increases with k s . The above inequalities (equations 18 and 
19) are satisfied for all k s . At early times, the data shows some violation of 
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the inequalities. This is an artifact of truncating the infinite Fourier series 
(Gibb’s phenomena) and is unphysical. At latter times when the influence 
of the high wavenumbers dies out, the inequalities axe satisfied by the scalar 
field. In any event it is clear that the range of the scalar concentration over 
which the pdf is non-zero decreases with time. 

In the presence of a turbulent field, the amplitude of various wavenum- 
bers decay in an unknown fashion. However, we still have the following 
inequalities, 

1 1 ±0 ° - 

^max(t) — ~ 5: n— ^ |<^(k, t)| (20) 

1 v 27T | fc |_! 

1 1 io ° ~ 

o' — $min(t) ^ ~/K = |<^(k, t)|- (21) 

1 v2tt | fc | =1 

The variance at time t is given by 

| ±oo 

»’(«)= ST £WM)f. (22) 

1*1=1 

Initially, <f> max and <f> m i„ are unity and zero, and the variance is some finite 
quantity. With mixing, the variance decreases, going to nearly zero at 
infinite time: 


1 ±oo 

cr 2 (t -> oo) = — £ [*(k, f = oo)) 2 - 0, 
2k l*f=i 


so that, 


As a result we have 


±oo 

|<^(k, t = oo) | -> 0. 

I*l=i 


— Oo) ♦ i 

. . .1 

ymin — OO) * ^ * 


(23) 


(24) 
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(25) 

(26) 



Hence, even in the presence of a turbulent velocity field we can expect the 
non-zero range of the pdf to decrease, contrary to the implications of the 
mapping closure solution of O’Brien and Jiang (1991). 

2.2 Conditional scalar dissipation 

As was seen in the Introduction, the evolution of the scalar pdf is determined 
by the conditional scalar dissipation, xM (equation 3). The mean scalar 
dissipation, e s , in terms of the conditional scalar dissipation, is 

- - S> = £ (27) 

While the conditional dissipation ratio, x(0)/ € *> P la y s an i m P ortant role 
in the evolution of the shape of the scalar pdf, the rate of the evolution is 

determined by £ 4 . 


Models. The models against which the HCS and DNS data are compared 
are: model Ml, the /3-pdf conditional dissipation ratio model of Girimaji 
(1992), and model M2, the model for normalized conditional dissipation, 
XW/X(0.5) (O’Brien and Jiang 1991). The model Ml is given by, 

XC0) _ » /*(! ~ Ap (28) 

e 4 <y 4 FW 


where /(?/>) is given by 


J(0) = / V '^{ln^-{ln<A)}+(l-M){ln( 1 -V’ , )-O n (l-^))}]W , )(V’-0 , )^ / . 

(29) 

The model M2 is given by 

-~~r = exp[— 2{er/ _ 1 (2i/> - l)} 2 ], (30) 

x(0-5) 
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where erf -1 represents the inverse error function. The most striking feature 
of the model M2 is that it is independent of the current scalar pdf. 

In Figure 4, the conditional dissipation ratios of HCS, DNS, and the con- 
ditional dissipation model of Girimaji (1992) are plotted against normalized 
mass fraction. The calculation of any high-order conditional statistic - like 
xi t s - from numerical data is usually fraught with statistical errors. Taking 
this into consideration, the agreement between HCS and DNS is good for 
the early stage ( a = 0.27) and adequate for the intermediate stage (cr = 0.1) 
of mixing. At the final stages of mixing, the conditional dissipation ratio 
of the HCS data is nearly constant over the range of mass fraction consid- 
ered. This is consistent with the fact that the HCS scalar pdf decays as a 
Gaussian at this range of variance. The DNS conditional dissipation ratio 
is minimum at tp = 0.5 and increases steeply for other values; this trend is 
inconsistent with the observed Gaussianity of the scalar pdf and is probably 
due to statistical error. The HCS data comes from a 128 3 simulation and, 
hence, can be expected to be statistically more accurate than the DNS data 
which comes from a 64 3 simulation. The model Ml agrees well with both 
HCS and DNS data during early and intermediate stages of mixing and 
shows the same trend as the HCS data at the latter stages. 

The above comparison between the model, HCS, and DNS data is only 
over a narrow range of mass fraction. In Figure 5, the normalized con- 
ditional dissipation x(V’)/x(0-5) of HCS is compared with the models Ml 
and M2 over the entire range of mass fraction. The models agree with HCS 
data during the early stages of mixing ( cr 2 > 0.05). During the latter stages 
(when <f> m ax < 1 and <j> m i n > 0), the agreement is not as good. The HCS 
values are set to zero outside the range (<^w„, (f> max )- In the middle of the 
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range, the agreement is adequate, whereas near the edges, it is poor. This 
disagreement near the edges of the range {(f> m i n ,<f>max) could be due to sta- 
tistical error in HCS data. Near the edges the probability of ^ is so low that 
at these values of mass fraction there are very few samples of scalar dissi- 
pation from which to compute the conditional dissipation. This statistical 
error makes it difficult to completely comprehend the evolution of condi- 
tional dissipation at the edges. However, it can be surmised from the earlier 
findings about <j> min (t) and <f> max {t) that the zero values of the conditional 
dissipation migrate inwards (from 0 and 1) with increased mixing. 

To circumvent the problem of too few statistical samples at the edges 
of the scalar range, in Figure 6, we plot the fractional scalar dissipation, 
*W)xM 0 f the HCS data and model Ml. The advantage of plotting this 
quantity is that when the statistical error is high due the low probability, the 
conditional dissipation is weighted down by the probability density itself. 
Here, the agreement between the HCS data and the ^-pdf-based model Ml 
is excellent. The jagged nature of the HCS-data plot at the earliest time 
is again due to the presence of very high wavenumber components in the 
scalar field. With time, the high wavenumber components die out quickly, 
yielding smooth statistics. 

Apart from the mathematical convenience it offers, the fractional scalar 
dissipation is a physically important quantity. It is this quantity, rather 
than the conditional dissipation alone, that determines the evolution of 
the scalar pdf (equation 3). In fact, the conditional dissipation always 
appears along with the probability density function (equation 27). This 
suggests that our efforts should be directed towards modeling fractional 
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scalar dissipation rather than conditional scalar dissipation. Moreover, the 
value of conditional dissipation outside the range ( <^ m ;„ , <f > max ) is irrelevant. 

Our inferences from Figures 4-6 are 

1. The conditional dissipation ratio, x(V0A»i calculated from the DNS, 
the HCS, and the model Ml agree adequately well. 

2. The normalized conditional dissipation, x(VO/x(0-5), of the models 
Ml and M2 agree well with HCS data during the early stages of mix- 
ing. At long times, the statistical error at the edges of the scalar range 
are too large to make a valid comparison. For this same reason, it 
is difficult to evaluate whether the normalized conditional dissipation 
always has a constant form (as suggested by the mapping-closure- 
based conditional-dissipation model of O’Brien and Jiang 1991), or 
if it tends to thin down with time, asymptoting to a spike shape at 
t{> = 0.5 (as suggested by Girimaji 1992). However, since <f) m i n and 
<f>max both tend towards the mean at long times (Figure 3), it appears 
that the latter is the more probable scenario. 

3. The fractional dissipation ratio, TbTxM ^ Q f JfCS and the model 
Ml show excellent agreement. Since it is this quantity that determines 
the scalar evolution, the agreement is of significance. 

2.3 Evolution of scalar variance and mean scalar dis- 
sipation 

The evolution of scalar variance (a 2 ) and the mean scalar dissipation (e s ) 
obtained from HCS for various initial length scales are plotted in Figure 
7. To precisely understand the effect of the velocity field on the evolution 
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of these quantities, we compare Figure 7 with the corresponding DNS data 
of Eswaran and Pope (1988, Figures 11 and 12). The biggest difference 
between the HCS and DNS data is in the evolution of e s . For all k s , the 
HCS-e s decreases monotonically in time, as expected. However, the DNS-e s 
shows a non-monotonic evolution for small k s . For k s < 4, the DNS mean 
scalar dissipation increases in the beginning. After peaking, the mean scalar 
dissipation decays, gradually at first and then more rapidly. This non- 
monotonic behavior is more pronounced for smaller k 3 . The implication of 
this on the variance evolution is that the DNS decay rate is higher than 
that of HCS for small k s . 

The above phenomenon, which is entirely due to the velocity field, is 
discussed in the next Section, where the effects of the velocity field are 
accounted for. 

2.4 Pdf of scalar gradient 

If the scalar evolves according to the heat-conduction equation, the scalar 
gradient also evolves according to the heat-conduction equation. However, 
the initial conditions of these fields are quite different. Due to discontinu- 
ities in the initial scalar field, the initial scalar-gradient field is composed 
of delta functions. In Figure 8, the evolution of the flatness-factor and 
superskewness of scalar-gradient component d(j>/dx\ is shown for k a = 8 
and k a = 16. For both cases, the moments approach Gaussian values at 
long times. The scalar-dissipation pdf (not shown) is, however, far from 
Gaussian. Eswaran and Pope (1988) find that the scalar dissipation is log- 
normally distributed. This log-normality is a consequence of the velocity 
field and is discussed in detail in the next Section. 
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3 Accounting for the effects of velocity field 


In this Section, we attempt to account for the effect of the turbulent velocity 
field on scalar mixing. Findings from the previous Section are invoked to 
effect simplifying assumptions in our analysis. 

The effect of the velocity field is most easily accounted for in the La- 
grangian frame of reference (x, t). Consider the evolution of the scalar 
field in a turbulent velocity field U(X, t). In this Section, we designate 
the Eulerian coordinates as X, rather than x, as was done in the previous 
Sections. Let the coordinate system evolve according to 

^h = U[X(x,i).(], (31) 


with the initial condition, 


x = X(0). 


(32) 


The velocity field U(X,t) evolves according to the Navier-Stokes equation. 
Now, let the Lagrangian coordinate system x be Cartesian. Then, the 
Eulerian system X is nonstationary, curvilinear, and nonorthogonal. In the 
coordinate system following fluid particles, the scalar concentration evolves 
according to 

dt dXidXi ’ ( j 

where D is scalar diffusivity. The molecular diffusion can be rewritten in 
terms of the Lagrangian scalar derivatives: 


d<f>(x, t) _ dx a dx b d 2 <f>(x, t ) 
dt ~ V dXidXi dx a dx b ' 
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The evolution equation of the pdf, G(^i,t), of the scalar concentration 
<f>(x,t) following fluid particles in an isotropic field can be derived using 
standard methods (Pope 1985): 

dG(^i, t) d 2 dx a dx b d(f>(x,t)d<f>(x,t) 

— (35) 

(In the absence of the velocity field, the above equation reduces to equation 
(3)). Due to the isotropy of the field, the Lagrangian scalar pdf Gi(ipi,t ) 
and the Eulerian pdf F(ip,t) are identical. 

Now we invoke the findings of the previous Section to simplify equation 
(35). It was found that, while the velocity field has a major effect on 
the time scale of scalar evolution, the pdf shapes during the evolution are 
unaltered for k a > 2 and only slightly modified for k s < 2. This, we argue, 
implies that the averaged effect of the velocity field conditioned upon the 
value of the scalar concentration is independent of the scalar value. That 
is, we suggest 

dx a dx b d<j>(x,t) d<f>(x,t) t dx a dx b d<t>(x,t)d<j)(x,t) 

( q xT—\<t> = i>l) « ( av flv )<— 51 1^ = Vlh 


K dXidXi dx a 


dx b 


K dXidXi ,x dx a dx b 


(36) 

where (f^fx 1- ) independent of the scalar field. This is not such a far- 
fetched assumption, since the value of the scalar concentration is a strong 
function of the initial scalar field, which is independent of the velocity field. 
Subject to the above assumption, equation (35) can be simplified as 

dG(ip t ,t) _ n/ dx a dx b ^ d 2 d<K x ,<) ^ = ^ 

dx a dx b 

(37) 


dt 


-D( 


dX t dx)'dind^ 
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Owing to the isotropy of the turbulence field and the initial scalar field, we 
can further simplify as follows: 


Then, 


dx a dx b 

dXi dXi 

dGjr/^t) 


&ab I dx c 9x c . 

T { dx,ox; ) = Sab{s) - 


(38) 


-(S)- 


d 2 




dt ' dtpidi^i 

where Xi( r fii) is the Lagrangian conditional scalar dissipation defined as 


(39) 




(40) 


dx a dx a 

Now we suggest that the effect of the velocity field on the evolution of 
(and, hence, on F(ift,t)) is only through the scalar ( S ), and that 
the Lagrangian conditional scalar dissipation is independent of the velocity 
field. Rescaling time as 


r = / \s(t))dt , (41) 

it is easily seen that the pdf evolution in the time coordinate r is inde- 
pendent of the velocity field. The velocity-field-independent Xi{i>i) i s the 
conditional scalar dissipation of a scalar field evolving under the influence 
of the heat conduction equation (1) alone, from the same initial conditions 
as the full problem. Thus, the effects of molecular action and velocity field 
can be decoupled. Clearly, these simplifications will be accurate for k 3 > 2, 
but less so otherwise. 


Mean scalar dissipation, e a . In isotropic turbulence, the conditional 
dissipation of a scalar field can be approximated as (from equations 36 and 
38) 


x(4>) = ( 


d<f> d<j> 

dXi dXi 


I <f> = ip) = (S)xi(^i)- 


(42) 
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The mean scalar dissipation is obtained by integrating equation (42): 


e ,(t) = (S(t))e,(t\ (43) 


where 


m n, W d * v 
£s{t) = 


K dX l dX % ' 

and e/(i) is the mean scalar dissipation in an equivalent heat-conduction 
system: 

^ &r a 5x a 

Some of the implications of this model are 

1. Case of zero diffusivity. In this case, there is no evolution of the 
Lagrangian field. Hence, remains constant. However, due to 

the velocity field, ( S{t )) will vary with time, and with it, 
will also change in time. 


2. Case of no velocity field. For this case, {S(t)) is always unity 
and, hence, the problem reduces to the heat conduction equation, as 
it should. 


3 ^ 0 and presence of velocity field. In this case, the turbulent 

scalar dissipation changes in time due to the effect of the velocity field 
( via change in (S(t))) as well as due to the decay of the Lagrangian 
scalar field ( via change in e/(t)). 

In order to understand the behavior of conditional scalar dissipation and 
mean scalar dissipation in turbulence, we need insight into the behavior of 
(S(t)) and e,(<). A study of the latter has already been performed in the 
previous Section. In the following subsection, we turn our attention to 

(S(t)). 
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3.1 Evolution of ( S(t )) 

The scalar S = 5 is a measure of the deformation of the coordinate 

frame caused by the velocity field. Differentiating equation (31) with respect 
to x a , we get 


dL _ d Ui(X(x,t),t) 

dt dx a 


dUi 

dXj 


(ja j 


(44) 


where 


_ ar,-(x,<) 

dx a 


(45) 


Equation (44) is the evolution equation of a material line element (Monin 
and Yaglom 1981). The quantity l ia stands for the i-th component of the 
material line that was initially coincident with the x a axis. Let I a represent 
the line-element vector whose components are l la , l 2a , and l 3a . Recall that 


dY,(x,Q) _ 
dx a 


(46) 


and, hence, I a (x, 0) = e a , the unit vector along the x a axis. 

The material line element l a is stretched (or shrunk) and reoriented by 
the velocity field. The line-element vector at any time t can be determined 
knowing the tensor B and the initial condition: 


1.(0 = B(0.e.. 


(47) 


The tensor B contains all of the one-point information of the fluid element 
and evolves according to (Monin and Yaglom 1981) 


dBa _ dU, 
dt dxk ky 


( 48 ) 
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What is S(t) in terms of / ta ? For evaluating S, we need the inverse of 
the matrix Let L, a be the inverse matrix; the components of this matrix 
in terms of Z to are 

Lu = -^7 [^22^33 — ^23^32] = y02 X = y 1 (^) 

1 1 ^12 
L \2 — ^[^32^13 — Z12Z33] = yfa X ^3)- e 2 = y ) 

1 1 C13 

L\3 = — [Z 12 Z 23 — ^22^13] = yfa X ^3)- e 3 = y 1 

1 1 C 2 I 

^21 = y [^23^31 — Z 33 Z 21 ] = y(^3 X ll)- e l = y > 

1 1 C22 

L 22 = y[hahl — Z 31 Z 13 ] = y(^3 X ll)- e 2 = y 1 

1 Is ^- x 23 

L 2 3 — y [hihz — W 23 ] = y (I 3 x li)-^3 = ~y ~ 5 

1 1 C31 

L 31 = "y [Z 32 Z 2 I — Z 31 Z 22 ] = y (U X ^2)- e l = y > 

1 1 C 32 

L32 — y[l 31^12 ~ ^32^1l] = y (1 1 X l2)*C2 = yT , 

1 1 C 33 

L33 — y[l 11^22 “ ^12^2l] = y Ol X l2)'63 = “y~ ’ 

where V is the ratio of volume (or density) of the fluid element. 



The numerator C xa can he identified as the i-th component of area of a 
surface element whose normal is initially coincident with the x a axis. The 
scalar S is then 


S(t) = ^Ci.(«)Ci.(t) = - ^\A](t) + A\(t) + AM, (51) 

where A a represents the magnitude of area of the material surface element 
whose normal is initially coincident with the x a axis. For an incompressible 
velocity field, the density ratio V is always unity. The compressibility of 
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the velocity field manifests itself on the scalar mixing through V. The 
implication, then, is that the compressibility only affects the time scale of 
scalar evolution, without modifying the shape of the scalar pdf’s. Indeed, 
Givi et al (1991) find that the 0 pdf replicates the scalar pdf in high-speed 
(compressible) mixing also. 

For an isotropic velocity field, the various area magnitudes (A\, A 2 and 
A 3 ) are all statistically equivalent. Hence, assuming isotropy, 

(5(()> = (A 2 (t)), (52) 

where A(t) represents the area at any time t of any typical material-surface 
element (of initial unit area) associated with the fluid element. 

A naive estimate of { S(t )) that ignores the effects of vorticity and strain- 
rate rotation can be made (a la Monin and Yaglom 1981). Subject to the 
above conditions (and incompressibility), we can expect 

A{t) ~ exp[(a x + a 2 )t], (53) 

so that 

(S(t)) ~ exp[2(a x + a 2 )f], (54) 

where a x and a 2 are the means of the positive and intermediate eigenvalues 
of the velocity gradient tensor. In real turbulence, however, the effects of 
vorticity and strain-rate rotation are important (Girimaji and Pope 1990). 

Girimaji and Pope (1990) have performed extensive studies of the evolu- 
tion of line, surface, and volume elements in statistically stationary, isotropic 
turbulence. One of their findings was that ln(A) has a Gaussian distribu- 
tion, which implies that A 2 has a log-normal distribution. Although no 
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direct estimate of { S(t )) was made, based on the findings of that study, we 
can surmise some aspects of the behavior of { S(t )). Let us define 

*> - ^ (55) 

so that we can write 

e s (t) ~ exp[j pdt]ei(t). (56) 

The exponential growth rate p(t) is a function only of the turbulence pa- 
rameters and is very likely to be positive at all times. In real turbulence, 
the temporal behavior of p(t) is likely to be very similar to the line and area 
logarithmic growth rates (Girimaji and Pope 1990). At time* = 0, p{t) will 
be nearly zero, owing to the lack of alignment between the material-element 
ellipsoid and the velocity-gradient tensor. As the alignment improves with 
time, p(t) would increase, peaking at about t m = 5r^ (r v , Kolmogorov 
time scale), when the alignment between the principle axes of the velocity- 
gradient tensor and that of the material-element ellipsoid is at the maxi- 
mum. Finally, as vorticity and the rotation of the strain-rate principal- axes 
reduce the alignment, p would asymptote to a lower value. Girimaji and 
Pope (1990) estimate the long-time growth rate of various moments of the 
area of surface elements and p(t) has the form 

( 57 ) 

dt T-q 

Implications of the behavior of { S(t )) on e s (<). The Lagrangian dis- 
sipation ei(t) is a function only of the initial length scale (/,) of the scalar 
field. Clearly, £/ decays monotonically in time, the rate of decay being larger 
for smaller l s and vice-versa. The non-monotonic behavior of e s in the DNS 
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data of Eswaran and Pope (1988) for large l a (small k 3 ) is, then, entirely 
due to ( S(t )). We can write 


din e a 

dt 


= P(t) + 


din e/ 
dt 


(58) 


If e a is to grow in time, we need p(t ) (which is always positive) to be greater 
in magnitude than the decay rate of In This can happen for small k a 
when the decay of e t is slow enough and when p(t) is near its maximum 
(t = t m ). At later times, when p{t) settles down at is smaller asymptotic 
value, € a will decrease due to the influence of ej. For the case of large k aj 
the decay rate of In q is always much larger in magnitude than p(t). Hence, 
for these values of k a , c a decreases monotonically. 


3.2 Pdf of scalar dissipation 


Eswaran and Pope (1988) find that at long times, the pdf of scalar dissi- 
pation is log-normal in the presence of isotropic turbulence. As mentioned 
in the previous Section, the pdf of scalar dissipation in the absence of the 
velocity field is not log-normal. We now investigate if our present analy- 
sis of the role of the velocity field can predict log-normality of the scalar 
dissipation. 

Let 0(= ■gxifx’) be the instantaneous value of the scalar dissipation. As 
before, we can write 

*(*) = T ab {t)e l ab {t), ( 59 ) 


where 




, rp _C±C ib 
dx a dx b ab V V ' 
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The matrix T is symmetric and positive definite. At long times, it would 
not be unreasonable to expect the relative orientation of the principal axes 
of T and the tensor 6 l ab to reach a stationary state. Then, 

e{t^oo)^\T(t)\\e\t)i ( 6 °) 

where \T\ and \0 l \ are the norms. This can be further simplified as 

e(t - oo) = ^\m, («) 

where A{t) again represents the area of a typical material surface associated 
with the fluid element. The physical meaning of this approximation is that 
at long times, the smallest length scale associated with a material element 
is the volume (of the fluid-element) divided by the area of a typical surface 
element. The scalar dissipation, being the square of the scalar gradient, 
should be inversely proportional to the square of the smallest length scale. 
So, at long times, 

ln(0) ~ ln(|0*|) + 21n(A) - 21n(F). (62) 

In incompressible turbulence, V(t) = 1 and, hence, ln(V) = 0. During the 
final stages of mixing, as the fluctuations in 10* | die out, the only fluctuations 
in the scalar dissipation are due to A(t). As mentioned before, in isotropic 
turbulence, the area of a surface element is log-normally distributed. Hence, 
the scalar dissipation at long times is also log-normal, and the log-normality 
is entirely the doing of the velocity field. 


29 



4 Mixing from non-standard initial condi- 
tions 

In Section 2 it was shown that, for the case of equal-quantity, non-premixed 
initial fields, the HCS scalar pdf and DNS scalar pdf are quite similar at 
all stages of decay. It is reasonable, then, to expect that the pdf’s of the 
heat conduction equation and full turbulence equation are similar for mixing 
from other initial conditions also, as long as the initial field is prescribed in a 
random manner. The effect of the velocity field on these non-standard mix- 
ing processes will be the same as that on the equal-quantity, non-premixed 
mixing. Recall that, throughout the analysis in Section 3, the only restric- 
tions on the initial scalar field are that the field is random and isotropic. 
The assumption in equation (36) is valid for any scalar field that meets the 
above requirements. 

In this Section, we study scalar fields evolving according to the heat 
conduction equation from non-standard initial conditions. Another impor- 
tant objective is to determine if the /9-pdf model replicates the pdf’s well 
for the case of non-equal, non-premixed initial fields. 

4.1 Non-equal, non-premixed initial fields 

In this subsection, we study the mixing of scalars from the following four 


initial conditions: 

Fi(ip, 0) = .85 + .155(1 - VO, H = 0.15, (63) 

F 2 (V>, 0) = M(1>) + .25(1 - VO, V = 0.20, (64) 

F 3 (iJ>, 0) - .75ty) + .35(1 - 0), n = 0.30, (65) 

F 4 (v>, o) = .65(^0 + .45(1 n = 0.40, (66) 
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with k s = 8. 

In Figure 9, the pdf of mass fraction at various stages of decay are 
plotted for the case of /x = 0.3. The HCS data is compared with the /3-pdf 
model at the corresponding value of mean and variance (Girimaji 1991a). 
The /3-pdf model appears to match the the data quite well over all stages 


of the mixing process. 

In Figure 10, the flatness-factor (S 4 ) of the mass fraction is plotted 
against variance for initial fields of various different fx s. Figure 1 1 represents 
a similar plot for the super-skewness (Se). With decreasing mean, the initial 
values of flatness-factor and super- skewness increase: 


S'n(O) = 


[m(i 


(67) 


For fx — 0.4 and ji — 0.3, S 4 and Se are monotonic functions of variance, go- 
ing to the Gaussian values of 3 and 15, respectively, at long times. However, 
for smaller /x’s (0.2 and 0.15), S 4 and S 6 increase with variance, overshooting 
the Gaussian values at the early stages and then decay slowly, asymptoting 
to their respective Gaussian values at long times. The smaller the mean, the 
more pronounced is the non-monotonic behavior. The /3-pdf-model values 
of S 4 and S 6 are also plotted in Figures 10 and 11. The model compares 
quite well with HCS data for S 4 and adequately well for S 6 . The model 
captures the non-monotonic behavior of the moments very well. 


Conditional scalar dissipation. The good agreement between HCS and 
the /3-pdf model begs the question, how does the conditional dissipations 
compare? In Figure 12, the fractional dissipation (FUVOxW/ 6 *) obtained 
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from HCS is compared with the model Ml (equation 28). The agreement 
is excellent over all stages of the mixing process. 

An important issue that has been brought up in the recent works of 
O’Brien and Jiang (1991) and Girimaji (1992) is whether the conditional 
scalar dissipation for non-equal, non-premixed mixing, is symmetric about 
ip = 0.5. The closure-mapping-based model of O’Brien and Jiang (1991) 
su SS es fs that the conditional dissipation is symmetric, whereas the model of 
Girimaji (1992) suggests that it is not. Although Girimaji (1992) presents 
qualitative reasons that the latter should be the case, no direct verification 
has been possible due to the lack of corresponding DNS data. In Figure 
13, the conditional dissipation calculated from HCS (for the case /i = 0.3) 
is plotted at various stages of scalar decay. As mentioned in Section 2, 
this statistic is more error-prone than fractional dissipation. Despite the 
unsmooth data, the following aspects of conditional dissipation are clear. 
At early stages of mixing, the conditional dissipation is symmetric about 
ip = 0.5. However, with progressing time, the conditional dissipation gets 
skewed, with the maximum migrating towards the mean (0.3 in this case), 
as suggested by the model of Girimaji (1992). The conditional-dissipation 
estimate in the vicinity of xp — 0.3 is quite reliable because of the high 
probability density in that neighborhood. As in the case of fi = 0.5, in this 
case also, we can expect <p min and <p max to migrate towards <p = fi with 
increased mixing. 

4.2 Partially-mixed initial field 

The main objectives of studying partially-mixed initial fields are (i) to de- 
termine if there is any qualitative difference in the long-time behavior and 
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(ii) to determine the degree of deviation from the /3-pdf model calculations. 
(The /3-pdf model cannot accommodate partially- mixed initial conditions.) 
The choice of partially-mixed initial fields is too numerous to cover all types 
of such fields. We consider the following three cases of initial pdf’s: 

F 1 (rP, 0) = 0.3<5(*/>) + 0.2 Sty - 0.25) + 0.2 Sty - 0.75) + 0.3^ - IP) 

F 2 {tl>,0) = 0.456(V>) + 0.21^ - 0.75) + 0.34£(t/> - 1), (69) 

F 3 0M) = 0.56(V-) + f7(0.5,l), (70) 

where U( 0.5, 1) represents the pdf of unit probability density uniform over 

the range (0.5, 1.). For these three cases, the means of the distributions are 

0.5, 0.5, and 0.35, respectively. Again, k a = 8 for all the simulations. 

In Figure 14, the pdf’s of case 2, F 2 (ip, t ), obtained from HCS are plotted 
at various stages of scalar decay. At early times, we have a triple-delta 
pdf. The delta-function at rj> = .75 is smeared out very rapidly at a rate 
disproportionate to the leveling out of the delta peaks at zero and unity. 
Very soon the scalar pdf appears like that of non-premixed initial fields (e.g., 
Figure 1). As expected, the /3-pdf model does not capture the intermediate 
peak in the early stages, but, as this peak is smeared out, the /3 pdf compares 
quite well with the HCS data. 

In Figure 15, the flatness-factors (£ 4 ) of the three cases are plotted 
against the variance. The model values of S 4 for the given value of mean 
and variance are also plotted in Figure 15. Figure 16 represents a similar 
plot of super-skewness. Our observations from Figure 15 and 16 are 

1. All the cases appear to tend to the Gaussian asymptotic state. 

2. The /3-pdf model agrees tolerably well with cases 1 and 2. 
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3. The agreement between case 3 and the corresponding /3-pdf model is 
not as good. For this case, the HCS moments S 4 and Se decrease 
with variance at the very early stages and then increase to Gaussian 
values. However, the model values grow monotonically to Gaussian 
values and are always higher than HCS values. 

5 Summary and conclusions 

This work represents an effort to improve our understanding of passive- 
scalar mixing in an isotropic turbulence velocity field. The most important 
contribution of this study is that it presents a simple, yet quantitative 
picture of turbulent mixing which adequately explains all of the known and 
observed behavior of such a mixing process. The findings of this study 
should simplify the modeling of turbulent mixing. The following is a brief 
summary of the findings. 

1. Many aspects of the scalar-pdf evolution in isotropic turbulence ob- 
served by Eswaran and Pope (1988) are the consequences of the effect 
of molecular diffusion on the initial random scalar field, rather than 
the effects of the velocity field itself. 

2. The effect of the velocity field is identified, isolated, and quantified. 
The velocity field manifests itself on molecular mixing via the quan- 
tity S(t) (equation 38) and only affects the time scale of the scalar 
pdf evolution. The analysis performed in Section 3 is valid for any 
velocity field that is isotropic and has a length scale larger than that 
of the scalar field. Hence, irrespective of whether the velocity field is 
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incompressible, compressible, or even evolving according to equations 
other than Navier-Stokes, we can expect the shape of the scalar pdf 
to evolve in the same universal manner. 

3. The following model for mean scalar dissipation is provided (rewriting 
equation 58) 

(n) 

where e is the mechanical dissipation of the turbulent field, v is the 
kinematic viscosity, and the second term on the right-hand side is 
a function of the initial spectrum and the diffusivity of the scalar 
field. This expression should provide a foundation on which realistic 
models for mean scalar dissipation can be constructed in the future. 
The model for e s that is currently in use is 

e, = CXp (72) 

where C\ is a numerical constant and k is the turbulent kinetic energy. 
This model is independent of the length scale or the diffusivity of the 
scalar field and, hence, unrealistic. 

4. An explanation for the observed log-normality of the scalar dissipa- 
tion is provided. A qualitative explanation for the log-normality was 
given previously by Hill and Bowhill (1978). Despite that, the present 
model represents an important contribution towards understanding 
turbulent mixing for, not only is it quantitative (relating the scalar 
dissipation to known turbulence parameters via material element de- 
formation statistics), but it is also a part of a larger model which 
incorporates all the observed properties of turbulent scalar mixing. 
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5. Analytical arguments and numerical evidence are presented question- 
ing the asymptotic correctness of the mapping-closure-based condi- 
tional scalar dissipation model of O’Brien and Jiang (1991), even for 
the case of molecular mixing (for which it is supposed to be most 
appropriate). 

6. It is demonstrated that even for the case of non-equal, non-premixed 
initial fields, the /?-pdf model performs quite well (Figure 9 - 11). Fur- 
ther, the /3-pdf-based conditional dissipation model (Girimaji 1992) 
agrees quite well with the HCS data (Figure 12). So, the conclusion is 
that the ft pdf is a good model for all non-premixed initial conditions. 

The methodology developed in this paper can be used to address several 
issues regarding chemical reactions in turbulent flows. Since chemical reac- 
tion is a process that is local in space, several of the processes that need to 
be modeled in turbulent combustion are also local in space and, hence, can 
be investigated in the absence of the velocity field. A study of a random 
scalar field of several species undergoing molecular diffusion and chemical 
reaction will be useful in assessing the effect of reaction on concentration 
fluctuations, the importance of temperature-concentration correlations, etc. 

Acknowledgements. I would like thank Dr. J. P. Drummond for his 
encouragement and support. This work was supported by the Theoretical 
Flow Physics Branch, Fluid Mechanics Division, NASA Langley Research 
Center, Hampton, VA 23665, under contract No. NASI - 18599. 


36 



References 


[1] Chen, H., Chen, S., and Kraichnan, R. H. (1989). Probability distribu- 
tion of a stochastically advccted scalar field. Physical Review Letters, 
63 (24), pp 2657 - 2660. 

[2] Eswaran, V., and Pope, S. B. (1988). Direct Numerical Simulations of 
the turbulent mixing of a passive scalar. Physics of Fluids, 31 (3), pp 
506 - 520. 

[3] Gao, F. (1991). An analytical solution for the scalar probability density 
function in homogeneous turbulence. Phys. Fluids A, 3 (4), 511 - 513. 

[4] Girimaji, S. S., and Pope, S. B. (1990). Material element deformation 
in isotropic turbulence J. Fluid Mech. 220 , pp 427 - 45b. 

[5] Girimaji, S. S. (1991a). Assumed [3 -pdf model for turbulent mixing: 
validation and extension to multiple scalar mixing. Combust. Sci. and 
Tech. 78 , 4 - 6, pp 169. 

[6] Girimaji, S. S. (1991b). A simple recipe for modeling reaction- rates 
in flows with turbulent combustion. AIAA-91-1792. AIAA 22nd Fluid 
Dynamics, Plasma Dynamics and Lasers Conference, June 24-26, 1991, 
Honolulu, Hawaii. 

[7] Girimaji, S. S. (1992). On the modeling of scalar diffusion in isotropic 
turbulence. To appear in Phys. Fluids. A. 

[8] Givi, P. (1989). Model-free simulations of turbulent reactive flows. 
Progress in Energy and Combustion Science, 15 (1), pp 1 ■ 10 i. 


37 


[9] Givi, P., Madnia, C. K., Steinberger, C. J., and Tsai, A. (1991) Large 
eddy simulations and direct numerical simulations of high speed turbu- 
lent reacting flows. State University of New York report, Department 
of Mechanical Engineering, SUNY, Buffalo, NY 14280. 

[10] Hill, R. J., and Bowhill, S. A. (1978). Lognormality of gradients of 
advected scalars. Physics of Fluids, 21 (6), pp 883 - 886. 

[11] Monin, A. S., and Yaglom, A. M. (1981). Statistical fluid mechanics } 
Vol. 2 (ed. J. L. Lumley), MIT Press. 

[12] O’Brien, E. E., and Jiang, T-L (1991). The conditional dissiption rate 
of an initially binary scalar. Physics of Fluids, 3, 3121 - 3123. 

[13] Pope, S. B. (1985). Pdf methods for turbulent reactive flows. Progress 
in Energy and Combustion Science, 11 , pp 119 - 192. 

[14] Pope, S. B. (1991). Mapping closures for turbulent mixing and reaction. 
Theoretical and Computational Fluid Dynamics, 2 , 255 * 270. 


38 



I 





Figure 1: Evolution of scalar pdf, F( t/>), for the case of /i = 0.5. Circles rep- 
resent DNS data, dashed line represents HCS data, and solid line represents 
/?- pdf model calculations. 










Figure 3: Temporal evolution of <f> m in and fima x for various initial length 
scales, (n - <f>mm{ks = l)i A - <j> max {k s = 2); Q — <f>min{ks = 2); O - 
= 2); + <f>„ linfl'i — 4); * - (f>max{^s = 4); X - <^ m in(^9 = 8); * 

- <f> m ax{ks = 8); solid line correspond to <f> ma x = 5 + ^e^p[ — i] and and the 
dashed line to <£ mtn = \ — |e:rp[-£].) 



Figure 4: Conditional scalar dissipation ratio vs. Normalized Mass Fraction 
for the case of p — 0.5. Comparison of DNS data, HCS data and model Ml 

(Girimaji 1991b). (<r = 0.27: — DNS, Q model, 0 HCS; cr = 0.10: 

DNS, □ model, ■ HCS; a - 0.02: DNS, A model, A HCS.) 
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Figure 5: Evolution of normalized scalar dissipation, xW/x(0-5), for the 
case of /t = 0.5. Circles represent model M2 (Gao 1991), dashed line repre- 
sents ITCS data, and solid line represents model Ml. 
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Figure 6: Evolution of fractional scalar dissipation, F(V>)x(V’)/ e *> for tlie 
case of /t = 0.5. Dashed line represents HCS data and solid line represents 
model Ml (Girimaji 1991b). 
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Figure 8: Evolution of IK'S flatness-factor and super-skewness of scalar- 
gradient d(j)/d. t,\ (fi = 0.5). (Flatness factor: O = 16; CD fc, = 8. 
Superskewness: A k., — 1(5; O k, = 8.) 
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1 igure 10. Scalar flatness-factor vs. scalar variance for non-'equal non 

premixed initial conditions, (ft = 0.15: □, HCS; /?- pdf model, ft = 

0.20: O (hexagon), HCS; /3-pdf model, ft = 0.30: O, HCS; — /?- pd 

model, ft = 0.40: O, HCS; £-pdf model.) 



Figure 11: Scalar super-skewness vs. scalar variance for non-equal non 
premixed initial conditions. (Symbols same as Figure 10.) 
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Figure 12: Evolution o 
case of fi = 0.3. Dashed 
/?-pdf model calculation 
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Figure 13'. Evolution of conditional scalar dissipation, for the case 

of /i = 0.3. HCS data. 
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Figure 16: Scalar super-skewness vs. scalar variance for partially-mixed 
initial conditions. (Symbols same as Figure 15.) 
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